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Abstract 

Recently several models and theories were proposed to describe and explain 
the scale-invariant behavior of crystal defect aggregations that possibly lead 
to self-affine geometry of fractured surfaces. Here the idealized general model 
of aggregate growth is proposed on the basis of the simple additive rules that 
correspond to one-step aggregation process. Some factors that can influence 
the self-affine nature of the aggregating system of solitary agents (monomers) 
and their aggregates (clusters) were brought to light: the rate of monomer 
particle exchange between clusters, cluster geometry, and initial cluster size. 
Two idealized cases were analytically investigated and simulated by Monte 
Carlo method in the Desktop Grid distributed computing environment to 
analyze "pile-up" and "wall" cluster distributions in different aggregation 
scenarios. The "pile-up" case with a minimum active surface (singularity) 
could imitate piling up aggregations of dislocations, and the case with a maxi- 
mum active surface could imitate arrangements of dislocations in walls. The 
change of scaling law (for pile-ups and walls) and availability of scale-free 
distributions (for walls) were analytically shown and confirmed by scaling 
analysis of simulated density and cumulative density distributions. The ini- 
tial "singular" symmetric distribution of pile-ups evolves by the "infinite" 
diffusive scaling law and later it is replaced by the other "semi-infinite" dif- 
fusive scaling law with asymmetric distribution of pile-ups. In contrast, the 
initial "singular" symmetric distributions of walls initially evolve by the dif- 
fusive scaling law and later it is replaced by the other linear scaling law with 
scale-free exponential distributions without distinctive peaks. The conclusion 
was made as to possible applications of such approach for scaling analysis of 
cumulative density distributions in experimental data. 
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1. Introduction 

In materials science the fractured surfaces demonstrate the self-affine ge- 
ometry on many scales. Moreover, the fractal analysis of fractured surfaces 
by projective covering and box-counting method shows that the fractured 
surface can be depicted not only by one fractal dimension, but also by mul- 
tifractal spectrum [H, 0). At the same time, surface roughness profiles of 
periodically deformed Al 0, 0, 0|, slip line morphology in Cu jg, 0], and 
dislocation patterns in Cu after tensile js|] also demonstrate the self-similar 
features on many scales. Recently, transition from homogeneous disloca- 
tion arrangement to scale-invariant structure was described by the statistical 
model of noise- induced transition ^ . Several other models and theories were 
proposed to explain the scale-invariant behavior of crystal defect aggregations 



that possibly lead to self-affine geometry of fractured surfaces [10 . 

The hierarchical defect substructures that were observed experimentally 
in deformed metals and alloys appear as a result of some aggregation pro- 
cesses among solitary crystal defects. It should be noted that many aggre- 
gation phenomena in other natural processes also take place by exchange 
of solitary agents (monomers) between their aggregates (clusters): phase 
ordering 11, 12 , atom deposition [l3j], stellar evolution [l3|, growth and dis- 
tribution of assets [lij], and city population even 16 1. In such aggregation 
processes monomers can leave one cluster and attach to another. Usually 
these exchange processes are described by an exchange rate kernel K(i, j), 
i.e. by the rate of transfer of monomers from a cluster of size i (detaching 
event) to a cluster of size j (attaching event). Generally, the rate of monomer 
exchange between two clusters depends on their active interface surfaces that 
are dependent on their sizes, morphology (line, plane, disk, sphere, fractal, 
etc), probability of detaching and attaching events, etc. 

Sometimes there is the preferable direction for exchanges, i.e. with asym- 
metric exchange kernels, K(i, j) ^ K{ j, i), like in coalescence processes 
in Lifshitz-Slyozov- Wagner theory [ill, IH, where big clusters "eat" smaller 



ones. The exchange rate kernel K(i, j) is defined by the product of the rate 
at which monomer detach from a cluster of size i and the rate at which this 
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monomer reach another cluster of size j. In Leyvraz-Redner scaling theory 
of aggregate growth [16j cities A of size i evolve according to the following 
rule: 



Ai + Aj K -H } A-i + A j+1 , (1) 

where K(i, j) is the exchange rate. That is, monomer (one person) leaves 
some of cities Ai of population i and arrive to some of cities Aj of population 
j. This can be considered as the generalized rule for the theory of growth 
and distribution of assets [15], if one can assume that Ai are persons with 
asset volume of i. 

Below the idealized general model of aggregate growth is proposed on the 
basis of this approach. The main aim of the work is to use the most profound 
features of aggregation kinetics and to find the simplest factors that can cause 
the observed self-affine properties of the aggregating system of solitary agents 
(monomers) and their aggregates (clusters). In this context, the numerous 
complex details of the real crystal defect aggregation processes will be hidden 
behind the idealized and simplified conditions only to emphasize the most 
general precursors of scale-invariant behavior of such complex systems. 

2. Model 

Here detaching and attaching processes are considered separately that in 
the general case could be characterized by different rates. Consequently, the 
different detach product kernel K d {n) = kdS (n) and attach product kernel 
K a {n) = k a S a (n) are taken into account, where k d and k a are the mea- 
sures of activation of detaching and attaching processes, n is the number 
of monomers in a cluster. In natural processes kd is usually determined by 
energy barrier for detachment from cluster and k a — by probability for at- 
tachment of migrating monomer to another cluster which in turn determined 
by kind of migration (instant hops from cluster to cluster, ballistic motion, 
random walking, or their combinations). Sd(n) = s d n a and S a (n) = s a n^ 
are the active surfaces of clusters, where a and j3 — exponents depending on 
the morphology of cluster (for example a = 1 for linear clusters and a = 2/3 
for spherical clusters, and a = (3 in the simplest case of clusters with the 
same morphology), s d and s a — the constants depending on the morphol- 
ogy of cluster and geometry of neighborhood (for example s d = 1 for linear 
aggregates and s d = v367r for spherical aggregates, and s d = s a = s in the 
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simplest case of clusters with the same morphology and neighborhood). The 
portion of clusters f(n,t) with n monomers at time t evolves according to 
the following equation: 



= K d (n + l)/(n + 1, t) + K a (n - l)f(n -l,t)- 
-K d (n)f(n,t)-K a (n)f(n,t) 

Actually this is the well-known master equation for so-called "one-step 
process", i.e. a continuous time stochastic Markov process, which range 
consists of integer n, and where transitions allowed only between adjacent 
integers |l7j]. 

The general model is based on the assumption that the migration time 
(movement of a monomer from one cluster to other) is much lower than the 
detachment (or attachment) time of a monomer from (to) a cluster. That is 
why there are no free monomers, and range of n is half-infinite (n = 2, 3, . . .). 
At the same time the total number of monomers is assumed to be constant 
and no generation sources exist. It should be taken into account that the more 
consistent system of evolution equations and more realistic coupling between 
migration and detachment (attachment) times should be used for rigorous 
comparison of this general model and each case of the aforementioned natural 
processes. 

In an asymptotic regime of high values of n equation (2) goes to: 

df (n,t) ^ d(Di (n)f(n,t)) + d 2 (D 2 (n) f (n,t)) 
dt dn dn 2 

which is the one-variable Fokker-Planck equation in general form with 
time-independent drift Di(n) = K d {n) — KJn) = n a s {k d — k a ) and diffusion 
D 2 (n) = K a {n) = n a sk a coefficients 0, The main difference between 
this formulation and other well-known models is the more general scenario 
with detaching and attaching as separate events is taken into account here. It 
allows us to express explicitly the availability of drift and diffusion terms, and, 
moreover, to bring to light the common symmetry properties and propose 



the ways to find the exact non-stationary solutions [20[ . 

The case a = corresponds to the clusters with the minimum active sur- 
face ("singularity"' of the active surface), namely for clusters with constant 
numbers of active monomers independent of the whole number of monomers 
n in it. For example, the ends of a line cluster can be its 2 active points 
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for detaching and attaching. Such configurations take place in queues (2 
active points), stacks (1 active point), etc. In the context of plastic deforma- 
tion of metals, "pile-up" aggregations of dislocations of a regular crystalline 
structure can be depicted by this scenario (see Fig. [T^,). 



In Out 
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Figure 1: Examples of (a) minimum active surface — pile- up of dislocations: it is possible 
to go in/out pile-up only through left/right ends of the pile-up; and (b) maximum active 
surface — wall of dislocations: it is possible to go in/out wall in any place. 

The case a = 1 corresponds to the clusters with the maximum active 
surface, namely for clusters where any monomer could be detached and free 
monomer could attached to any place. Again, in plastically deformed met- 
als, "wall" aggregations of dislocations of a regular crystalline structure can 
be depicted by this scenario (see Fig. [lb). In Leyvraz-Redner scaling the- 
ory of aggregate growth of city population [l^| it can correspond to linear 
dependence of arrival or departure rates as a function of city population. 

The case < a < 1 corresponds to the clusters with the bulk volume 
of monomers shielded by active surface. For example, circumference of disk, 
perimeter of fractal, outer layer of sphere can be the correspondent active 
surfaces for detaching and attaching. In solid state physics, such configura- 
tions take place in regular arrangements (compact like voids and spare like 
fractals) of point-like defects of a crystalline structure. 

The case kd = k a corresponds to the equiprobable activation of detaching 
and attaching processes, for example, for periodic tensile deformation it can 
correspond to the equiprobable activation of detaching and attaching dislo- 
cations from opposite sides of wall. Finally, it means the absence of the drift 
term in (3). 

2.1. Pile-up - minimum active surface 

For the case a = clusters have the minimum active surface and under 
condition kd = k a one can get the well-known heat equation, where a = sk a , 
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df(n,t) _ d 2 f(n,t) 

dt ~ a On* ' [) 

which has numerous particular solutions that are dependent on the initial 
and boundary conditions. In the reality kd ^ k a , and the aggregation process 
is described by the homogeneous heat equation with space-independent drift 
and diffusion coefficients. It is well-known fact that it leads to "a diffusive- like 
kinetic universality class". 

Below for illustrative purpose we consider the time evolution of the ini- 
tial singular distribution of clusters of the same size (n ) that exchange 
the monomers. It is actually the first boundary value problem for domain 
< n < oo, / (0, t) — 0, / (n, 0) = S (n — n ), where 5 (n) is a Dirac delta 
function. This problem has the well-known solution: 



f(n,t) 



2\/nat 



exp 



[n - n ) 
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exp 



K n + n ) 
Aat 



(5) 



In practice, / (n, t) could be experimentally determined for high values 
of n and t. That is why / (n, t) is of interest for large values of n ^> n and 
t > n 2 /4a: 



f(n,t) = 



nriQ 



■ exp 



n 

4ot 



(6) 



2^(atf 2 

The essential point is that on the initial stage of this one-step aggre- 
gation process the probability distribution function / (n, t) does not "feel" 
the boundary condition / (0, t) = 0. From the physical point of view the 
cluster size cannot reach the 0-boundary (and equation cannot "feel" the 
boundary condition) immediately, but only after some time to, when the first 
cluster disappear due to detaching monomers. It means that on the initial 
stage t < t the aforementioned problem formulated as the classic Cauchy 
problem for infinite domain — oo < n < oo with the same initial condition 
/ (n, 0) = 5 (n — n ) and the well-known solution: 



/t«o {n, t) 



2^/nat 



exp 



(n - n Q y 



Aat 



(7) 



and for n ^> n it goes to the following shape: 
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(8) 



Thus we have two different scaling laws for the probability distribution 
functions for pile-ups / (n, t) for high values of n ^> n : 



for t < to ( "infinite" diffusive scaling law) 

f t<t0 {n,\t)^\- 1/2 f{\- 1/2 n, t), 

for t > to ( "semi-infinite" diffusive scaling law) 

f t>t0 (7i,At)->A-7(A- 1/2 ra, t). 



(9) 



(10) 



From practical point of view it is not easy to distinguish scaling laws 
(Q and ([TO]) and difference between them by analysis of probability density 
functions f(n) for various moments of time (Fig. [2^), especially from scarce 
and noisy experimental data. 
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Figure 2: Pile-up size distributions: probability density functions f(n) (a) and their cu- 
mulative distribution functions CDF(n) (b) for different stages of aggregation kinetics (j4|). 
Thickness of lines corresponds to the increasing values of t in the order from the thinnest 
to thickest line: (a) 100; 1000; 2000; 4000; and (b) 10; 10 2 ; 10 3 ; 10 4 ; 10 5 ; 10 6 ; 10 7 . 

But it is much more easier to observe the scaling laws by analysis of the 
cumulative distribution functions. For example, the visual comprehension of 
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scaling law fllUp can be easily obtained by viewing the steady shift of the 
correspondent cumulative distribution functions (CDFs) to the right side in 
Fig. [2b. The shift becomes diffusive (i.e. proportional to y/i) after t > 
10 4 , when the first cluster will reach the 0-boundary, i.e. by continuous 
detachments some cluster will decrease its size to and disappear. 

It means that one can determine the qualitatively different scaling regimes 
- "infinite" and "semi-infinite" ffTUl) — in this aggregation model with 
additive rules by the scaling analysis of the probability distribution function 
/ (n,t) for high values of n ^> uq. As it will be shown below in section HJ 
scaling analysis of cumulative density functions is more efficient for prac- 
tical purposes. The unbounded case (Cauchy problem for infinite domain 
— oo < n < oo) corresponds to a situation with the constant number of 
clusters (monomers only redistribute between them), and the bounded one 
(first boundary value problem for domain < n < oo) — to a situation with 
irreversible decrease of the number of clusters. 

2.2. Wall - maximum active surface 

For the case a = 1 clusters have the minimum active surface and under 
condition kd = k a one can get: 

df (n, t) d 2 (nf(n,t)) 

~^ = a M ' (U) 

which has numerous particular solutions / (n, t) ~ F (n/t) that are depen- 
dent on the initial and boundary conditions. Again for illustrative purpose 
we consider the time evolution of the initial singular distribution of clusters of 
the same size (no) that exchange the monomers. It is actually the first bound- 
ary value problem for domain < n < oo, f (0, t) = 0, / (n, 0) = 5 (n — n ) 
with the following solution: 



ft A V"° 
/ ( n > *) = -T-7= ex P 
atyfn 



(n + np) 
at 



/2^/nraA 



(12) 



Again, on the initial stage of this one-step aggregation process the prob- 
ability distribution function / (n, t) does not "feel" the boundary condition 
/ (0,t) = 0, because from the physical point of view the cluster size cannot 
reach the boundary immediately, but only after some time to, when the first 
cluster disappear due to detaching monomers. When t to and n ~ n , 
the solution of ( fTTj) has not influenced by differences in n significantly and 
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that is why it has approximate solution, which is close to solution of the heat 
equation: 




(13) 



But for the later stage t > t solution of ( ITT]) will be close to: 



ft>to {n, t) -> -*—= 




(14) 



and for n ^> uq it goes to the following shape: 



ft>t (n, t) -> 




a 



~ n ' 

exp - 

L at. 



(15) 



Thus, again we have two different scaling laws for the probability distri- 
bution function of walls / (n, t) for n 3> n : 

• for t < to ("infinite" diffusive scaling law) 



Again, as in the case with pile-ups, it is not easy to distinguish scaling 
laws fTTHT) and (fTTj) and difference between them by analysis of probability 
density functions f(n) for various moments of time (Fig. [3^l), especially from 
experimental data. 

And it is much more easier to observe the scaling laws by analysis of 
the cumulative distribution functions (Fig. [3b). The visual comprehension 
of scaling law (fT7|) can be easily obtained by viewing the steady shift of the 
cumulative distribution functions (CDFs) to the right side in Fig. [3b. The 
shift becomes linear (i.e. proportional to t) after t > 10 2 (but not after 
t > 10 4 , like it was for pile-ups!), when the first cluster will reach the 0- 
boundary, i.e. when some cluster by continuous detachments will decrease 
its size to and disappear. But the most essential point is that shortly after 
t > t the size distribution for wall configuration becomes exponential (|T5|) . 
that is scale-free one and without any "apparent" peak value. 



f t<t0 {n,\t)^\- 1/2 f{\- ll2 n, t), 



(16) 



• for t > to and n ^> n ( "semi- infinite" linear scaling law) 



ft>to (n, At) -»■ A 



/(An, At). 



(17) 
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Figure 3: Wall size distributions for different stages of aggregation kinetics (a) and their 
cumulative distribution functions (b). Thickness of lines corresponds to the values of t in 
the order from the thickest to thinnest line: (a) 1; 10; 20; 40; 80; 160; and (b) 1; 10; 10 2 ; 
10 3 ; 10 4 ; 10 5 



One of drawbacks of the simplified model is that the equations 01]) and 
( ITT]) and their solutions (j5J) and ffl2]) give the idealized and rough represen- 
tation of the aggregation kinetics. Tin fact, solutions (jSJ) and ffl2]) allows 
the "long-tails" of distribution, i.e. non-zero density of clusters for arbitrary 
high or low values of n in the ranges n > n Q and n < n Q for any stage of 
evolution. For example, any initial cluster of size n = 10 2 will need at least 
t = 10 2 detachment steps to disappear. But the analytical solution (5) gives 
the non-zero values of cluster distribution f(n,t) for any small n, i.e. does 
not prohibit immediate disappearance of this initial cluster even at the first 
step after start of the aggregation process. In a real physical process it is 
impossible due to physical limitations of diffusive character of the one-step 
aggregation process among clusters with pile-up morphology. This situation 
is well-known and actively investigated in the class of so-called "first-passage 



problems" 17]. That is why, the analytical solutions and predictions of 
the model cannot be directly applied to the correspondent practical situa- 
tions without taking into account the limits of their physical counterparts. 
Here, we investigated the aforementioned scaling laws and transitions be- 
tween them by the following simulations. 
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3. Simulation 



The two aforementioned primitive cases of cluster aggregation were simu- 
lated by Monte Carlo method to illustrate the different cluster distributions in 
different aggregation kinetics. The numerous initial configurations of clusters 
with various numbers of monomers in each of them were used in simulations 
with a conserved number of 10 6 aggregating monomers as a whole. Below the 
results on initial configurations of 10 4 clusters with 10 2 monomers are consid- 



ered (and the other vast configurations will be reported later 20(). Further 
the analytical calculations in section [2] will be compared with the results of 
simulation and for this purpose the number of MCSs will be assumed to 
equal to the number of time steps t. All simulation runs were performed in 
the distributed computing infrastructure (DCI) "SLinCA@Home" (Scaling 
Law in Cluster Aggregation) (http://dg.imp.kiev.ua/slinca) on the basis of 



BOINC SZTAKI Desktop Grid (DG) technology |2l|, |22 



3.1. Pile-up — minimum active surface 

The "pile-up" aggregation of monomers (see Fig. [T^,) has a minimum 
active surface and was simulated under condition that kd ^ k a . The kinetics 
of rearrangements from the initial configuration is shown in FigJH 

After some Monte Carlo sweeps (MCSs) of the initial configuration the 
broad peak appears, which corresponds to the average cluster size. The size 
distributions of pile-ups (Fig. @Ji,b) become asymmetric, but they preserve 
their one-peak shapes. And their CDFs (Fig. Hfc,d) demonstrate the "diffu- 
sive" scaling law, which becomes evident after some number of MCSs, namely 
after t > 10 4 . It can be seen visually by the shift of the simulated CDFs to 
the right side (proportionally to y/i) on the log-log plot in Fig. EH and can be 
compared with the same shift of the analytically calculated CDFs in Fig. [3b. 
The two straight lines are given in Fig. Hb as a guide for eyes to emphasize 
that f(n, t) ~ n for n < n Q and t > £ , that is follows from (5). 

3.2. Wall — maximum active surface 

The "wall" aggregation of monomers (see Fig. [lb) has a maximum active 
surface and was simulated under the same condition that kd ^ k a . The ki- 
netics of rearrangements from initial configurations is shown in Fig. [5j It can 
be seen that size distributions of walls (Fig. |5^,b) become asymmetric also. 
Moreover, they do not preserve their one-peak shapes after some number of 
MCSs and become scale-free, as it was predicted by the analytical solution 
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c) d) 

Figure 4: Size distributions of pile-ups for different stages of aggregation kinetics in: (a) 
linear and (b) double logarithmic coordinates (the straight lines are given as a guide for 
eyes to emphasize linear part of f(n,t)); and their cumulated distribution functions for: 
(c) linear and (d) logarithmic values of cluster size. Figures near symbols denote the 
number of Monte Carlo steps (MCSs). 
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(TPoj) for t > t . After some Monte Carlo sweeps (MCSs) the initial peak, 
that corresponds to the initial cluster size n , disappear and the scale-free 
distribution without distinctive peak appears. This transition is much more 
evident in the log-log representation of cluster size distributions on Fig. [5b. 
Their CDFs (Fig. [5b, d) also demonstrate the "linear" scaling law (IT7|) . which 
becomes evident after some number of MCSs, when the shift of CDFs to 
the right side of the plot in Fig. [5H becomes linearly proportional to t after 
t > 10 2 . One can easily compare it with the similar shift of CDFs in Fig. |3b). 
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Figure 5: Size distributions of walls for different stages of aggregation kinetics in: (a) 
linear and (b) double logarithmic coordinates; and their cumulated distribution functions 
for: (c) linear and (d) logarithmic values of cluster size. Figures near symbols denote the 
number of MCSs. 
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4. Discussion of Results 

Comparison of analytical and simulation results shows that monomer ag- 
gregation (at least, in pile- up and wall configurations of clusters) can be 
satisfactorily described by the proposed theoretical model of one-step ag- 
gregation processes. The important thing is that some scaling laws can be 
determined experimentally by scaling analysis of size distributions (and cor- 
respondent probability density functions — PDFs) and CDFs of the real 
aggregating ensembles. Below some attempts of such analysis are shown on 
the basis of aforementioned simulation results. 

4-1. Pile-up — minimum active surface 

The application of "infinite" diffusive scaling law <Q to simulated size dis- 
tributions (Fig. |6k) and CDFs of pile-ups (Fig. [6b ) shows the clear difference 
between two scaling regimes: "infinite" OH]) and "semi-infinite" (JTUJ) : 




a) b) 

Figure 6: Different stages of aggregation kinetics of pile-ups represented by scaled: (a) 
size distributions and (b) cumulated distribution functions. Figures near symbols denote 
the number of MCSs. 



• The "infinite" diffusive scaling law (Q corresponds to the initial stage 
t < to, i.e. the scaled size distribution curves are attracted to the initial 
symmetric size distribution; and the scaled CDF curves — to the initial 
straight line on probability plot in Fig. [6b (straight line on such plot 
corresponds to CDF of normally distributed values). 



14 



• The "semi-infinite" diffusive scaling law ( II Op takes place later (for t > 
t ) after some transient period, i.e. the scaled size distribution curves 
are attracted to some other asymmetric size distribution; and the scaled 
CDF curves — to other curved line on probability plot in Fig. [6b. 

The approximate transition region between two scaling laws can be ob- 
served in Fig. as the solitary "water line" curve for t = 10 4 (denoted by 
down-directed triangles (v) i n Fig- and by left-directed triangles (<) in 
Fig. |6b), i.e. for the moment of the first-passage of a cluster through the 
point of disappearance (n = 1). 

4-2. Wall — maximum active surface 

The application of "infinite" diffusive scaling law (I16p to simulated size 
distributions (Fig. [Tfet) and CDFs of walls (Fig. |8^) shows the other clear 
difference between two scaling regimes: "infinite" diffusive f lT6|) and "semi- 
infinite" linear ( ITT)) : 




Cluster Size Cluster Size 



a) b) 

Figure 7: Wall size distributions for different stages of aggregation kinetics scaled by: (a) 
"infinite" diffusive scaling law ([Tt?| : (b) "semi- infinite" linear scaling law (|17p . Figures 
near symbols denote the number of MCSs. 



• The "infinite" diffusive scaling law (fT6)) corresponds to the initial stage 
t < to, i.e. the scaled size distribution curves are attracted to the initial 
symmetric size distribution Fig. [TJi; and the scaled CDF curves — to 
the initial straight line on probability plot in Fig. 
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• The "semi-infinite" linear scaling law (ITTj) takes place later (for t > t ) 
after some transient period, i.e. the scaled size distribution curves 
are attracted to some other asymmetric and scale-free size distribution 
(on Fig. [7b); and the scaled CDF curves — to other curved line on 
probability plot (Fig. |8b). 




a) b) 

Figure 8: Scaled cumulated distribution functions for walls scaled by: (a) "infinite" diffu- 
sive scaling law (fT6|) ; (b) mixed scaling laws — "infinite" diffusive (p~6|) and "semi- infinite" 
linear (|17j) (see marks in legends of markers) . Figures near symbols denote the number of 
MCSs. 

The transition between two scaling laws can be clearly observed in Fig. [8b, 
where two family of collapsed curves are seen. 

Scaling analysis of simulated (or experimentally obtained) size distribu- 
tions (and PDFs) is very hard and prone to errors task, because of the some 
uncertainty related to proper and optimal binning procedures for preparation 
of histograms. In contrast, the scaling analysis of experimental and simu- 
lated CDFs has much higher stability to statistical deviations and outliers 
among data. The better collapse of scaled CDF curves could be seen even by 
visual comparison of scaled simulated histograms (Fig. Eh and Fig. [TJi) with 
their scaled CDFs (Fig. [6b and Fig. [8b). 

The initial "singular" configurations of walls after some time sweeps 
evolve to scale-free distributions (Fig. Uh), that is different from pile- up dis- 
tributions with broad peaks for the same time (Fig. [6^). That is why in 
practical sense the notion of "average cluster size" seems to be meaning- 
less for this kind of cluster size distributions (i.e. for exponential-like wall 
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arrangements) without distinctive peaks. This observations can give some 
insights as to possible roots for self-affine arrangements of defects and their 
manifestations at the surface of plastically deformed metals. Because in the 
asymptotic limit of high values of n > n and t ^> t Q the exponential distribu- 
tion of wall sizes f|T5|) can be approximated as f(n,t) ~ A(t) — B(t)n, where 
are A(t) and B(t) are some functions of time steps. That is why in practice 
to estimate the level of hierarchy of defect substructures and related scale 
of damage it is necessary to measure not only the integral "averaged" char- 
acteristics of defect ensembles (density per area, average size, etc.), but also 
to define more specific characteristics, related with their size distribution, 
like CDFs and their parameters j24|. One of the examples of such scaling 
analysis applied for CDFs was demonstrated for statistical analysis of fea- 
tures observed on surfaces of plastically deformed Al single crystals under 
real-time video monitoring and processing (2o1 |. 

At the moment some details of the simulation process cannot be compared 
directly with the theoretical results for large n and t because simulations were 
carried out for relatively low number of aggregating monomers (~10 6 ) and in 
the limited time range (< 10 4 steps). That is why to check the exact solutions 
and make more reliable conclusions the bigger simulations are under way now 
in the DG DCI "SLinCA@Home" . 

In the more general sense, scaling analysis of experimental distributions 
could allow us to determine their intrinsic symmetry and bring to light the 
corresponding aggregation scenario for many other applications related with 
aggregation phenomena, for example for investigation of city population dy- 
namics 26 



5. Conclusions 

The idealized general model of one-step processes to characterize aggre- 
gate growth was proposed that allowed us to find the scaling characteristics 
of some aggregation scenarios. Some factors that can cause scaling transition 
and appearance of the self-affine size distribution of the aggregating system 
of solitary agents (monomers) and their aggregates (clusters) were brought 
to light: the rate of monomer exchange between clusters, cluster geometry, 
and initial cluster size. It was shown that the simplified aggregate growth 
can be described by the one-variable Fokker-Planck equation in general form 
with time-independent drift and diffusion coefficients. In some aggregation 
scenarios it could be transformed to the equivalent equations, from which 
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the non-stationary analytical solutions with some initial and boundary con- 
ditions can be obtained. 

Two primitive cases of cluster aggregation were considered analytically 
and simulated by Monte Carlo method to illustrate the different cluster distri- 
butions in different aggregation kinetics: model with minimum active surface 
(singularity) that corresponds to "pile-up" aggregation of dislocations, and 
model with maximum active surface that corresponds to "wall" aggregation 
of dislocations. It was shown that initial "singular" symmetric distribution 
of pile-ups evolves according to "infinite" diffusive scaling law and later it is 
replaced by the other "semi-infinite" diffusive scaling law with asymmetric 
distribution of pile-ups. In contrast, the initial "singular" symmetric dis- 
tributions of walls initially evolve according to the diffusive scaling law and 
later it is replaced by the other linear scaling law with scale-free exponential 
distributions without distinctive peaks. 

From the practical point of view and in the more general sense, the scaling 
analysis of experimental distributions of aggregating monomers and clusters 
will allow to determine their intrinsic symmetry, bring to light the corre- 
sponding aggregation scenario, and influence of aforementioned parameters 
(rate of monomer exchange between clusters, cluster geometry, and initial 
cluster size). 
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